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Abstract 

We have studied the detection of long-lived staus at the IceCube neutrino 
telescope, after their production inside the Earth through the inelastic scatter- 
ed ■ ing of high energy neutrinos. The theoretical predictions for the stau flux are 
j — ■ calculated in two scenarios in which the presence of long-lived staus is naturally 
associated to viable supersymmetric dark matter. Namely, we consider the cases 
t— i ■ with superWIMP (gravitino or axino) and neutralino dark matter (along the 
CN ' coannihilation region). In both scenarios the maximum value of the stau flux 
qq ' turns out to be about 1 event/yr in regions with a light stau. This is consistent 
with light gravitinos, with masses constrained by an upper limit which ranges 
# ^ ' from 0.2 to 15 GeV, depending on the stau mass. Likewise, it is compatible with 
axinos with a mass of about 1 GeV and a very low reheating temperature of order 
100 GeV. In the case of the neutralino dark matter this favours regions with a 
low value of tan (3, for which the neutralino-stau coannihilation region occurs for 
smaller values of the stau mass. Finally, we study the case of a general super- 
gravity theory and show how for specific choices of non-universal soft parameters 
the predicted stau flux can increase moderately. 



1 Introduction 



The existence of high energy neutrinos has been proposed in several theoretical mod- 
els. They are produced in astrophysical sources via collision of hadrons or of hadrons 
with the surrounding photons. Since neutrinos are not deflected by magnetic fields and 
do not lose energy through interactions with the background while travelling from the 
sources to us, they carry all the relevant information about the nature of the astrophys- 
ical sources. Detection of such neutrinos has been attempted in current experiments 
like AMANDA and will also be pursued in future kilometer scale detectors such as 
IceCube at the south pole and KM3NeT at the Mediterranean Sea. 

Through inelastic scattering with nucleons (N) inside the Earth, these high energy 
neutrinos can produce exotic particles which, if charged and long-lived, may be de- 
tected in the above mentioned km 3 Cerenkov detectors, thereby opening a window to 
new physics. This idea was first proposed in [1] and further studied in [21 0] within 
the context of a supersymmetric theory. In particular, they considered the process 
u + N — > q + l, assuming a spectrum where the lighter stau, Ti, is a long-lived NLSP to 
which the squark and slepton promptly decay. In addition, the possibility of producing 
and detecting these long-lived particles within the Universal Extra Dimension model 
has been recently explored [3] . In [H [3] the number of staus that might be detected in 
IceCube was estimated by a Monte-Carlo simulation for fixed slepton, stau and wino 
masses and for three different squark masses. In [21 H] an analytic estimation of this 
number was performed for the SPS7 supersymmetry benchmark point and for a toy 
model with sparticle masses right above the experimental limits. An interesting re- 
lated possibility, namely the production of long-lived staus after the collision of high 
energy cosmic rays with nuclei in the upper atmosphere, was studied in [6] showing that 
staus arriving at large zenith angles could be detected in IceCube. Finally, stau detec- 
tion taking into account other neutrino sources which act as background for neutrino 
telescopes has also been explored [7j. 

In this work we have further explored the possibility of producing staus inside the 
Earth and detecting them at neutrino telescopes. More specifically, we have probed the 
parameter space of the Minimal Supersymmetric Standard Model (MSSM) from the 
point of view of a supergravity theory where parameters are defined at the GUT scale. 
Using the renormalization group equations to calculate the resulting supersymmetric 
spectrum, we calculate for each point the theoretical predictions for the stau flux at 
IceCube. Furthermore, we have also investigated the implications which arise when 
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a supersymmetric solution to the problem of dark matter is also imposed. Indeed, 
the lightest supersymmetric particle (LSP), when neutral, is an excellent dark matter 
candidate in R-parity conserving models. We have contemplated two scenarios which 
provide an LSP dark matter candidate as well as a long lived stau when it is the 
next-to-lightest supersymmetric particle (NLSP). In the first of them, the LSP is the 
lightest neutralino [Sj. In order to have a long lived enough NLSP in this tight 
degeneracy between its mass and that of the stau is necessary. The second scenario 
is that with a super WIMP LSP, namely the gravitino [9] or the axino [10] • Both 
particles are characterized by extremely weak interactions, which entails a long lived 
NLSP. In exploring these two scenarios, we have taken into account the most recent 
experimental constraints (such as limits on the masses of supersymmetric particles, 
and on low energy observables) , together with the present bounds on the relic density 
of cold dark matter. 

This paper is organized as follows. In Section 2 we describe the physical processes 
leading to slepton detection at neutrino telescopes. In Section 3, we describe in detail 
the two above mentioned scenarios, namely that where the neutralino is the LSP and 
that with a superWIMP LSP. We then show the resulting theoretical predictions for 
stau detection rate in both of them. We also explore the possible enhancement of 
the stau flux in scenarios with non-universal soft supersymmetry-breaking parameters. 
Our conclusions are left for Section 4. 

2 Slepton detection at neutrino telescopes 

The existence of ultra high energy cosmic rays as well as the detection of TeV photon 
emissions from galactic and extragalactic sources (supernovae remnants, active galaxy 
nuclei and gamma ray bursts) are a strong indication for neutrino emission from the 
same sources. Waxman and Bahcall (WB) estimated an upper bound for this flux 
assuming a proton cosmic ray flux proportional to E~ 2 , motivated by first order Fermi 
acceleration models [11]. Mannheim, Protheroe and Rachen (MPR) also determined 
an upper limit for diffuse neutrino sources [12] in almost the same way as WB, but 
instead of assuming a specific cosmic ray spectrum, they defined their spectrum based 
on current data at each energy. The resulting neutrino flux is approximately one order 
of magnitude above the WB prediction. Notice that the bounds on the neutrino flux 
could soon become more constrained by IceCube data [IB"] . 
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High energy neutrinos reach the Earth at a given point on its surface, As 
they propagate through the Earth towards the detector, and due to Standard Model 
interactions, their flux will be attenuated. Denoting the flux at the Earth's surface as 
F®(E U ), then the flux at a position I reads 

F,(£„) =F 9 {E„) e* V { f p{l\e) aSM{K) dl'X , (2.1) 

Ui 9 m P ) 

where p(l, 9) is the Earth's density at each point, m p is the proton mass, and a 
corresponds to the neutrino-nucleon scattering cross-section calculated in the Standard 
Model P3]. At this point, the neutrino interacts with a nucleon, resulting in the 
production of a pair of supersymmetric particles, which will further decay to staus and 
cross the remaining distance to the detector, placed at I = 0. 

To calculate the number of events at the detector per unit time and unit area, 
we have to multiply the parton level supersymmetry cross-section (a SUSY ), the corre- 
sponding parton distribution function (PDF), f(x,Q 2 ), the density at each point, and 
the neutrino flux taking into account its attenuation. Upon integration in I, cos6, x, 
Q 2 and neutrino energy, the resulting flux (usually expressed in yr _1 km -2 ) reads 

r 1 r l(B oil 9) r 1 r Q ™* r E ™ ax da SUSY 

N = 2n dcosO / f^-i-dl / dx dQ 2 / dEJ(x,Q 2 ) F t (E v ) (2.2) 

Jo Jo m p Jx min JQ 2 min JEi nin dxdQ 

The integration limits in x, Q 2 are given by kinematics, as well as the lower limit of 
neutrino energy. For an upper limit we have chosen E™ ax = 10 10 GeV. Taking higher 
limits will not affect the results, given the E~ 2 dependence of the flux. The limits on 
I depend on three points: 

• Energy losses: The average energy loss of a particle traversing a column depth z 
(where dz = p(l, 6)dl) is given by 

-(%)=a + 0E, (2.3) 

where E is the energy of the particle, a describes ionization energy losses and (3 is 
the radiative energy loss, which receives contributions from bremsstrahlung, pair 
production and photonuclear scattering. The parameter a is nearly constant as 
a function of the mass of the particle, a~2x 10~ 3 GeV cm 2 /g . 

For low energies (E -C E cr = a//3), the range of the stau is dominated by either 
ionization energy loss or by the lifetime of the particle, and it scales linearly with 
energy. 
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On the other hand, (3 shows a dependence with both the energy and the mass 
of the particle. In particular, there is a 1/m^ dependence with stau mass, and 
an increase of /3 with stau energy. In this work we have used a parametrization 
of P obtained in [15] by a Monte-Carlo evaluation of the stau range including 
electromagnetic energy losses^. 

The dependence of (3 with the mass of the particle is also true for muons, and 
the relation /3 M m M « /^m^ holds. This means that (3^ is approximately three 
orders of magnitude bigger than f3?, causing the muon range to be much smaller 
than the stau range. This crucial fact, allows stau detectability even though 
their production cross-section is about three orders of magnitude smaller than 
standard model processes p]. 

• Stau lifetime: If the stau is not long-lived enough, it may decay before it reaches 
the detector. We will further comment on this issue when we discuss in more 
detail the two scenarios we have studied. 

• Track separation: The angle between the two staus has to be wide enough to 
ensure the detector to be able to discriminate the arrival of the two particles 
from a one-particle event. We therefore require that the particles are at least 
50 metres apart at arrival. On the other hand, if the angle is too wide, one of 
them, or both, may miss the detector. Hence, we demand them to be at most 1 
km apart from each other. 

In principle, muons produced by upgoing atmospheric neutrinos could mimic the 
signal of individual staus traversing the detector. This kind of background is elliminated 
with the requirement that two simultaneous tracks have to be observed in the detector. 

This leaves the production and subsequent detection of a muon pair, as 
the main source of background. The stau track separation is a key point in order to 
discriminate stau pair flux from this signal [3] . Due to their shorter range, the detected 

1 The effect of weak interactions on the parametrization of 0, especially those coming from charged 
current interactions can be comparable to that of electromagnetic interactions [16] . The impact 
on the range is maximal for pure left-handed mass eigenstates. It increases with the energy, and 
decreases with the mass of the stau. Given the E~ 2 shape of the neutrino flux, detected staus are 
expected to come from neutrino interactions close to the threshold energy for squark production 
(i.e., approximately I0 6 GeV). At that energy, these weak interaction effects are not important, and 
hence, in [16j it is concluded that event rate estimates with energy losses parametrized as in [15] are 
reasonably reliable. Moreover, the areas of the MSSM parameter space that we have explored give 
rise to a lighter stau with a small left-handed component. 
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muons are only those produced in the vicinity of the detector. Hence, the separation of 
the tracks of a pair is smaller or of the order of 100 metres. On the contrary, staus 

can be produced at much larger distances and as a consequence their track separation 
can be as large as several hundreds or even thousand of metres. Thus, although the 
dimuon flux largely exceed the stau flux, for track separations above 100 m there should 
not be any significant contribution from this background. Moreover, making use of the 
energy deposition of the events it is possible to further reduce the dimuon background. 
All these issues are explained in detail in Ref. [3]. 

Hence, the limits on the range z(cos9, I) are calculated as follows. The upper limit, 
Zmax, corresponds to the minimum between the particle's range due to energy losses 
and the distance it can travel before decaying, I = 7CT. The lower limit is determined 
by the requirements on the track separation: for a given angle, I has to be such that 
A = I tan 9 lab is between 50 m and 1 km. 

The analysis of particle tracks in neutrino telescopes is calibrated for muons. The 
energy of the incoming muons that reach the IceCube detector can be reconstructed 
through their energy loss AE per column depth Az, provided that > E C J '. According 
to equation ( 12. 3ft . AE / Az = a + /S^E^. Thus staus would be detected as muons with 
reduced energy Ed = E? m^/m?, since /3f = (3^ m^/nif. We will use this reduced energy 
to calculate the total rate of events at IceCube, as the acceptance of the telescope highly 
depends on detected energy, as well as on the incoming direction of the particles. Note 
that the critical energy for staus is much higher than that for muons and so for arrival 
energies below E~ r ~ 10 5 GeV, one has AE/ Az m a and hence it would not be possible 
to estimate their energies. 

Finally, in our calculation we have used the effective area for upward going muons, 
determined by the IceCube collaboration [I?] , averaged in the angular direction as it 
is detailed in Table 1 of Ref. [2J. This implies a suppression factor in Eq. (12.21) . which 
is smaller than 1.25 for detected energies above 1000 GeV and which can be as large 
as a factor 10 for energies below 100 GeV. We have computed the effect of weighting 
with the IceCube effective area for a representative set points in the parameter space 
and have found it to induce a reduction in the total flux of approximately a factor of 
1.5. This can be qualitatively understood from the fact that the detected stau energy 
distribution peaks at approximately 1000 GeV [2J. 
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3 Results 



We are now ready to determine the theoretical predictions for the flux of staus that 
could be observed in neutrino telescopes. In doing so, we will work within the context 
of a supergravity theory, in which the set of soft supersymmetry-breaking parameters 
are considered as inputs at a high energy scale, which we will take as the scale at which 
gauge couplings unify, Mqut ~ 2 x 10 16 GeV. The renormalization group equations 
(RGEs) for the MSSM are then numerically solved to evaluate these parameters at the 
electroweak scale where the supersymmetric spectrum is calculated. The minimization 
of the Higgs potential leaves the following condition to be satisfied by the Higgsino 
mass parameter at the SUSY scale, 

„ —mjj tan 2 (3 + m 2 Hi 1 

a = Hu 9 ' ^ - -Ml, 3.4 

P tan 2 /3-l 2 Zl K J 

in terms of the ratio of the vacuum expectation values of the Higgs doublet, tan/9 = 
(H®)/(H%), which is also considered an input in our analysis. Notice that the sign of 
fj, is also left undetermined. 

Inelastic scattering on nucleons is the dominant interaction process of high energy 
cosmic neutrinos both in the atmosphere and the inside of the Earth. The products of 
this interaction are sleptons and squarks which can ultimately decay into the lighter 
stau, when this is the NLSP and long-lived enough. The leading processes in the MSSM 
involve exchange of charginos and neutralinos along a t-channel. The corresponding 
expressions of the differential cross-section for neutrino-quark scattering into a pair of 
sleptons and squarks are explicitly shown in the Appendix. 

The theoretical predictions for the resulting stau flux are obviously dependent on the 
initial structure of the soft terms, which is a function of the (yet unknown) mechanism of 
supersymmetry breaking. In the following study we will start by assuming universality 
of the soft-parameters at the GUT scale, studying the so called Constrained MSSM 
(CMSSM). The CMSSM is fully specified by a common scalar mass, m, a common 
gaugino mass, M, a trilinear parameter, A, the sign of the \i term and tan/3. 

In exploring the supersymmetric parameter space we will impose several experi- 
mental constraint in order to guarantee phenomenological consistency. More specif- 
ically, we will consider the LEP bounds on the masses of supersymmetric particles, 
as well as on the lightest Higgs boson. Moreover, we will also include the current 
experimental bound on the branching ratio of the b — > 57 decay, which sets the 
most stringent constraints in the scenarios we analyse. In particular, we will impose 
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2.85 x 10~ 4 < BR(6 — > S7) < 4.25 x 10 -4 , obtained from the experimental world average 
reported by the Heavy Flavour Averaging Group [13], and the theoretical calculation 
in the Standard Model |19j, with errors combined in quadrature. 

As already explained in the introduction, we will consider two main scenarios in 
which the presence of long-lived staus is well motivated and associated to solutions to 
the problem of the dark matter. In particular, we will start by analysing the CMSSM 
scenario with a supersymmetric superWIMP (gravitino or axino LSP) and then we 
will extend our study to the case of neutralino dark matter in the coannihilation re- 
gion with the stau. Finally, we will explore the effect of non-universalities in the soft 
supersymmetry-breaking parameters on the predictions for the stau flux. 

3.1 superWIMPs 

Two possible superWIMPs [2111 121] are viable dark matter candidates within the frame- 
work of supersymmetric theories, namely the gravitino and the axino. In both cases 
the stau, when it is the NLSP, is long-lived (easily exceeding = 1CT 9 s) due to the 
smallness of the couplings (gravitational and Peccei-Quinn scale, respectively) govern- 
ing its decay into the LSP. Therefore, both situations are suitable frameworks that 
would give rise to a flux of staus in neutrino Cerenkov detectors. Let us first briefly 
comment on both possibilities. 

On the one hand, the gravitino, when it is the LSP in a supergravity scenario, can 
be an excellent candidate for dark matter [20] • Gravitinos can be thermally produced 
during the reheating of the Universe, by ordinary processes involving scattering and 
decays of particles in the primordial plasma. Besides, a non-thermal production class 
of process also exist when the NLSP has a lifetime such that, being shorter than 
the age of the Universe, it is still long enough so that it decouples from the plasma 
before decaying to the gravitino LSP. Which contribution to the relic density is more 
significant depends on the reheating temperature and on the gravitino mass, but in 
general both have to be considered [221 [22]. 

In late decays of the NLSP into the LSP, additional electromagnetic and hadronic 
showers are produced. If the decay takes place after Big Bang Nucleosynthesis (BBN), 
the products of these showers may alter the abundances of light elements [21]. More- 
over, late injection of electromagnetic energy may distort the frequency dependence 
of the cosmic microwave background spectrum from its observed blackbody shape 
[251 1261 [27] . Preventing these effects leads to constraints on the supersymmetric pa- 
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rameter space that have to be imposed in addition to the usual experimental bounds. 
Recently, it has been pointed out [2H] that metastable charged particles can form bound 
states with light nuclei, opening new channels for thermal reactions and enabling the 
catalyzed BBN. In this case 6 Li and 9 Be overproduction becomes an issue and places 
very strong upper bounds on staus lifetimes, ry < 5 x 10 3 s (see, e.g., [291 130]). 

These constraints exclude extensive areas of the parameter space. In particular, in 
the case of the CMSSM and for moderate gravitino masses, they disfavour the regions 
with neutralino NLSP [3TI |32| [33] , leaving only some regions in which the stau is the 
NLSP. Therefore, the presence of long-lived staus in the case of gravitino LSP is very 
natural. Indeed, in this scenario, the stau decays to the gravitino and a r lepton at 
tree level, via gravitational interactions with a lifetime [3"H 151]. 

On the other hand, the axino (fermionic supersymmetric partner of the axion) is 
another possible dark matter candidate [2"T] . Due to the smallness of its coupling to 
ordinary matter, 1/ f a (with /„ ~ 10 11 GeV being the Peccei-Quinn scale), it has similar 
properties as those of the gravitino. For example, the axino relic density also receives 
contributions both from thermal and non-thermal production processes. An important 
difference is, however, that the lifetime of the NLSP is considerably smaller (since the 
axion coupling is much smaller than the gravitational one). Hence, the NLSP typically 
decays before BBN and the parameter space is generally free from constraints on the 
abundance of light elements. 

We proceed now to analyse the detection properties of the stau in neutrino tele- 
scopes when it is the NLSP in the case of either the gravitino or axino dark matter sce- 
nario in the CMSSM framework. We will therefore analyse those points in the CMSSM 
where the stau is the lightest observable supersymmetric particle (LOSP) assuming that 
the rest of the parameters (gravitino or axino mass and reheatingtemperature) can be 
chosen in order to achieve viable gravitino or axino dark mattery. 

Thus, for each point in the parameter space we have calculated, using the expres- 
sions in Appendix |A], the sfermion production from neutrino inelastic scattering inside 
the Earth. We have then used the code ISAJET [35J to check the subsequent decay 



2 This might not be possible for all the points in the parameter space, but given a specific scenario 
for gravitino or axino dark matter the corresponding BBN constraints and regions with correct DM 
relic density can easily be superimposed on our plots. 
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Figure 1: Supersymmetric spectrum corresponding to the CMSSM with M = 1000 GeV, 
m = 20 GeV, A = and tan/3 = 10. 

chains and branching ratios of these parent sleptons and squarks into the lighter stau 
and we have made an estimation of the average energy of the parent particle carried 
by the staus. In particular, we have found that the following relations 

w \ (Ej) ; (Ef) « i {Eg) , (3.6) 

hold for the entire region of the CMSSM parameter space where the stau is the LOSP 
with M < 2 TeV. This is similar to previous estimations performed for the SPS7 bench- 
mark point [2]. For completeness, we show in Fig.[T]a representative spectrum obtained 
in the CMSSM for M = 1000 GeV, m = 20 GeV, A = and tan/3 = 10. 

Our choice for PDF's in Eq. (12.21) corresponds to those obtained from the Cteq6PDF 
package [36] . However, in order to evaluate the dependence of the resulting flux on the 
particular choice of PDF's, we have repeated our calculations using also the MRST2004 
package [HZ]- We find a very small variation in the predicted stau flux, that we can 
quantify as approximately a 2%. Moreover, in the computation with the Cteq6PDF 
package, we have also evaluated the uncertainty in the resulting stau flux which is due 
to the PDF's error. Following the procedure detailed in Ref.|36J we have found that 
the error is approximately a 2%. 
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Figure 2: Theoretical predictions for the stau flux at IceCube as a function of the stau mass 
for different values of tan in the CMSSM scenario with A = 0. 



Using then the procedure sketched in the previous Section, we have computed 
the flux of staus at the IceCube detector. We have explored a slice of the CMSSM 
parameter space, setting A = 0, fi > 0, and varying the scalar and gaugino mass 
parameters, m and M, in the range — 2 TeV, retaining only those points where the 
stau is the LOSP and the various experimental constraints are fulfilled. 

In Fig. [21 we represent the stau flux as a function of stau mass for both WB 
and MPR neutrino fluxes and for tan/? = 10 and 50. Let us first examine the case 
tan (3 = 10. As expected, the production of heavier staus is suppressed, and therefore 
the flux decreases with the stau mass. We obtain, however, a band and not a line 
when we plot flux versus mass. This can be explained in the following way. Neutrino 
interactions leading to stau pair production involve the production of a squark, whose 
mass determines the energy threshold of the whole process. Squark masses are strongly 
dominated by gluino masses, which are roughly proportional to the common gaugino 
masses. Therefore, for the same value of the stau mass, the points on the (m, M) 
plane corresponding to lower values of M give rise to a lighter spectrum and thus lower 
energy thresholds. This, together with the E~ 2 shape of the neutrino flux, implies a 
higher number of stau pairs for lower values of M. Turning to the figure, it is now 
easy to understand that, given a value for m?, the highest value of the flux corresponds 
to the lowest value of M for which the stau is the LOSP. The lowest value, in turn, 
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corresponds to the highest value of M, i.e., that corresponding to m = 0. In the figure, 
the lower limit on the stau mass corresponds to the LEP experimental constraint on 
the Higgs (stau) mass for tan/5 = 10 (50). 

For the tan j3 = 50 case, the flux is lower than that for tan /3 = 10 in approximately 
one order of magnitude. Notice that when tan/3 increases, the off-diagonal elements of 
the sparticle mass matrices become larger and the L-R mixing becomes more important. 
This implies a larger mass gap between the two mass eigenstates of the third family. 
Thus, for a fixed stau mass, the rest of the spectrum is heavier for tan (3 = 50 than for 
tan/3 = 10. This, according to the previous threshold energy argument, implies a lower 
flux. Notice finally that for tan/3 = 50 the experimental constraint on the branching 
ratio of the b — > s'j, which implies a lower bound on the common gaugino mass, leads 
to an upper bound on the predicted stau flux of about 0.2 event yr _1 for the MPR flux. 

Thus the most optimistic results are obtained for light staus, and the lowest possible 
stau mass in the tan (5 = 10 scenario sets an upper limit to the predicted flux of about 
1 event yr" 1 when the MPR flux is used (with the WB prediction this is reduced to 
0.3 event yr _1 ). Notice that these results are smaller than those obtained in [3] due 
to the more constrained framework of supergravity analyses. Despite the fact that, 
as explained in the previous Section, the stau events can be distinguished from the 
dimuon background, these results evidence that, at best, several years of data from 
IceCube would be necessary to claim a positive signal. We will later see how this is 
modified when non-universal soft parameter are included. 

In order to understand which regions of the CMSSM parameter space could be 
explored using this technique, in Fig. [3] we superimpose the theoretical predictions for 
the stau flux on the (m, M) plane for A = and tan (3 = 10, 50. The most optimistic 
predictions regarding detectability occur for low values of M, corresponding to the 
regions with lighter staus and the flux decreases as the gaugino mass increases, with 
almost no dependence on the scalar mass parameter (when the stau mass is RG-evolved, 
the main contributions to it come from terms depending on gaugino masses). When 
the MPR prediction for the neutrino flux is used, a significant increase of the resulting 
stau flux is obtained. 

To sort out the limits on the stau lifetime imposed by catalyzed BBN, the mass of 
the gravitino is constrained by an upper limit, 



which can be deduced from Eq. (13.51) . setting the stau life equal to 5 x 10 3 s. It can then 




(3.7) 
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Figure 3: (m, M) plane of the CMSSM scenario with A = and tan (3 = 10, 50. The region 
below the thick solid line corresponds to that with stau LOSP. Coloured regions represent 
the theoretical predictions for the stau flux at IceCube in the case of WB neutrino flux, 
whereas vertical solid lines correspond to the results when the MPR flux is used. Gridded 
areas are experimentally excluded by LEP constraints on the Higgs and stau masses, whereas 
the ruled area is excluded by b — > 57. The white area on the lower left corner of the plot 
for tan (3 = 50 is excluded due to the occurrence of negative stau squared-mass. 



be easily seen that the maximum gravitino mass ranges approximately from 0.2 GeV to 
20 GeV, depending on the mass of the stau. For such light gravitinos, the relic density 
is fully dominated by thermal production, 



nTP h ^ Q27 ( T R \ / 100GeV \ frrvM \ 



(3.1 



where rrig(fi) is the running gluino mass [38, 39J. A relic density of about 0.1, in 
agreement with WMAP data [3D], can be recovered by an appropriate choice of the 
reheating temperature, depending on the mass of both the gravitino and the gluino. 
As Tr oc mg/m|, we can make a rough estimation of the maximum value of Tr as a 
function of the stau mass, 

r ^ 5 ' 2xl ° 7 (lo^) I/2GeV - < 3 ' 9 > 

Similar upper bounds for the reheating temperature have been derived [23, [2HI I3TJ and 
used to study its possible determination at the LHC [31]. It is finally worth pointing out 
that in this case the contribution to Qdm from non-thermal production is negligible. 
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As in the case of the gravitino, axino thermal production is a function of the re- 
heating temperature [2D 112] • As emphasized in [13] , the correct relic density can be 
obtained for nearly any point in the (M, m) by fixing the axino mass and the reheating 
temperature. In particular, the regions with a small value of the gaugino mass (where 
the stau flux is maximal) could correspond to an axino with mass m^ ~ 0(1 GeV) 
which must have a small reheating temperature, Tr ~ 200 GeV, in order to satisfy the 
constraint on its relic abundance. Notice that although a gauge invariant treatment of 
the axino thermal production [12] is only valid for larger Tr, this value can be obtained 
within the context of a quitessential kination scenario [Sj. BBN constraints in this 
case are easily fulfilled since the stau lifetime is usually smaller than 1 second. 

3.2 Neutralino LSP 

Let us now consider a second interesting class of MSSM scenarios, those with a neu- 
tralino LSP as the dark matter candidate. Under these circumstances, in order to have 
a sufficiently long lifetime, the stau NLSP needs to be almost degenerate with the light- 
est neutralino. Interestingly, the quasi-degeneracy of the stau and lightest neutralino 
is also welcome from the point of view of neutralino dark matter. In such a situa- 
tion the relic abundance of neutralinos is largely suppressed through a coannihilation 
mechanism [13] , which makes it possible to find agreement with the constraints on the 
dark matter abundance [46]. Thus, this is another example in which supersymmetric 
scenarios which solve the dark matter problem can naturally provide long-lived staus. 

The dependence of the stau lifetime on Am = m^ — m^o and on the left-right 
content of the mass eigenstate was studied in [47J. It was shown that in the regions 
where Am > m T , the 2-body decay r — > t\° is allowed, and it is the dominant process. 
In this case, the stau decays very rapidly, with a lifetime smaller or of the order of 
10~ 17 s, and therefore never reaches the detector. On the other hand, if Am < m T , this 
channel closes, and typical lifetimes become (9(10 _6 s) or greater. This is long enough 
for staus to reach the detector before decaying 

This is indeed a very restrictive requirement which only leaves a very narrow al- 
lowed band in the parameter space. In fact, such a small mass- difference between the 
neutralino LSP and the stau implies a too small relic density for the neutralino when 
its mass is small. The presence of staus with a sufficiently large lifetime is therefore 

3 For example, given the mean free path for a stau, = ct? E^/m?, imposing I? > 1 km for a stau 
of ~ 150 GeV and an energy of 10 6 GeV, this would imply a lifetime of at least 10~ 9 s. 
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Figure 4: (Am, m~) plane of the CMSSM scenario with A = 0, tan/3 = 10,50 and 
neutralino LSP. Coloured regions represent the theoretical predictions for the stau flux at 
IceCube in the case of WB neutrino flux, whereas solid vertical lines correspond to the 
results when then MPR flux is used. Regions where WMAP relic density is reproduced are 
represented by gridded areas. 

only compatible with neutralino dark matter above a certain mass scale. 

In order to illustrate this, we have represented in Fig. H] the predicted flux of staus 
at the IceCube detector in the (Am, m?) plane for an example of the CMSSM with 
A = 0, fi > and tan/3 = 10, 50, and for both a WB and a MPR neutrino flux. 
As commented above, the flux vanishes for Am > 2 GeV. For Am < 2 GeV the flux 
increases as the stau mass decreases and can be as large as 2 event yr -1 for the MPR 
flux when the stau mass is close to its experimental lower bound. The WB flux is 
represented by means of a colour code, whereas the regions corresponding to different 
values of stau flux assuming a MPR neutrino flux are delimited by solid lines. The 
regions compatible with neutralino dark matter are represented on the same plane by 
means of gridded areas (for those points where the WMAP relic density is reproduced). 

Compatibility of neutralino dark matter with observable staus is only possible for 
tan j3 = 10 and for m^ « 300 GeV, thus implying a stau flux between 0.1 and 0.5 events 
yr _1 , in the optimistic MPR case or between 0.01 and 0.1 event yr _1 if we work with 
a WB neutrino flux. For larger values of tan/3 the region where the neutralino relic 
density is compatible with WMAP results is shifted towards larger values of the stau 
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mass and therefore is associated to much lower values of the stau flux. Once more, in 
spite of the good background discrimination, these small fluxes would require, at best, 
several years of data from IceCube. 

3.3 General supergravity with non-universal soft terms 

In the former Sections we have performed our calculations within the framework of the 
CMSSM. We will now briefly explore how the theoretical predictions for the stau flux 
vary when a more general supergravity scenario with non-universal soft parameters is 
considered. More specifically, we will try to identify possible non-universal schemes 
that lead to an increase of the theoretical predictions for the stau flux. 

As we discussed above, the production rate of staus is very sensitive to the value of 
the squark masses, the resulting flux increasing in the presence of light squarks [U [3j. 
Likewise, a decrease in the slepton masses is also welcome in order to increase their 
production cross-section. Moreover, as we can see in the Appendix, the expressions 
( 1A.12I) corresponding to the chargino-mediated sfermion production are proportional 
to the factor \Z^ Zl fc | 2 . This implies that these contributions are enhanced when the 
lighter chargino is a pure wino state (i.e., \Z}^\ ~ 1). 

A decrease of the low energy values for the squark masses, relative to the (right- 
handed) stau mass can be obtained if the ratio of the gluino and bino mass parameters, 
M 3 /Mi, decreases at the GUT scale. The decrease in the squark masses induces a 
reduction of the fi parameter through the increase of the positive contributions to 
the RGE of m 2 H u - This in turn implies an unwanted enhancement of the mixing in 
the chargino mass matrix. Thus, in order to compensate for the decrease in the /z 
parameter and have a pure wino as the lighter chargino, as well as reducing its mass, 
we must also decrease the M 2 parameter at the GUT scale. Smaller values of M 2 also 
imply a decrease in the masses for left-handed sleptons. Although, as mentioned above, 
this is also potentially good to increase their production cross-section, one should bear 
in mind that eventually the sneutrinos can become almost degenerate with the lighter 
stau and be long-lived. In that case left-handed squarks would generally cascade down 
to the lightest neutralino and this, in turn, mostly to sneutrinos which could propagate 
through the Earth without decaying into staus, thereby significantly attenuating the 
stau flux. We have avoided this situation by making sure sneutrinos decay promptly 
(with a lifetime smaller than 10 -9 s), thus constraining the non-universality in M 2 . 

Finally, we can also attempt to increase the fi parameter, thus further enhancing 
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Figure 5: Supersymmetric spectrum corresponding to non-universal soft mass parameters 
as in Eq.( l3~T0l with M = 1000 GeV, m = 20 GeV, A = and tan/5 = 10. 

the wino composition of the lightest chargino. This can be done by introducing a 
departure from universality in the Higgs mass parameters at the GUT scale. More 
specifically, we have considered an increase of the value of w/jj with respect to the 
mass of the rest of the scalars. 

As an specific example, following the above arguments we have taken the following 
set of non-universalities in the gaugino and scalar masses, 

M 2 = 0.6 Mi 
M 3 = 0.5 Mi 

2 2 
m H D ,Q,U,D,L,E = 171 ) 

m 2 Hu = 0.5 m 2 . (3.10) 

The resulting supersymmetric spectrum is depicted in Fig.[5]for M = 1000 GeV and 
m = 20 GeV, clearly evidencing the decrease in the masses of the squarks and heavy 
sleptons relative to the stau mass with respect to the universal case discussed in the 
previous Sections and shown in Fig.[U 

Given the change in the SUSY spectrum, we have redone the analysis of the decay 
chains that produce staus from squark and sleptons. We have checked that in this 
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example the relation between the resulting stau energy and the parent squark and 
slepton reads 

(E 7 ) « 0.8 (Ej) ; (Ek) « 0.5 (E ? ) . (3.11) 

The increase in both quantities with respect to the universal case (see Eq. ( 13.61) ) is due 
to the smaller mass differences among the supersymmetric particles which cascade to 
the lighter stau, which leads to a smaller suppression according to Eq. (IA.14l) . These 
relations approximately hold throughout the whole region of the parameter space with 
stau LOSP. 

The resulting theoretical predictions for the stau flux are represented in Fig. [U] as a 
function of the stau masses for both the WB and MPR flux for A = and tan (3 = 10. A 
moderate increase is observed with respect to the universal case of Fig. El For example, 
the results using the WB flux can reach now values of almost 1 event yr _1 , while those 
corresponding to the MPR flux are approximately a factor three larger. 



4 Conclusions 



We have explored the possibility of detecting long-lived staus at neutrino telescopes, 
after their production inside the Earth through the inelastic scattering of high energy 
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neutrinos. More specifically, we have calculated the production rate of a pair of staus 
in terms of the high energy neutrino flux, evaluated their energy losses after traversing 
the distance to the experiment, and taken into account the track separation of the stau 
pair and the detector efficiency in order to compute the theoretical predictions for the 
resulting stau flux. We have studied two generic scenarios in which long-lived staus 
are naturally associated to a supersymmetric solution to the problem of dark matter. 

On the one hand, we have considered the case of superWIMP dark matter in which 
the LSP is either the gravitino or the axino. Exploring the areas of the CMSSM in 
which the stau is the lightest observable particle we observed that the number of stau 
pairs is bounded to be below 1 event yr _1 when the MPR neutrino flux is used. These 
predictions decrease by approximately a factor four for the WB flux. The largest values 
of the stau flux correspond to small values of tan /3 and to regions with a small value of 
the common gaugino mass. These areas of the parameter space can be consistent with 
viable gravitino dark matter for light gravitinos (mg < 1 GeV), in order to avoid the 
stringent BBN constraints, or axinos with mass ~ (9(1 GeV) and a small reheating 
temperature, Tr ~ 200 GeV. 

On the other hand, we have explored the case with a neutralino LSP. In this sce- 
nario, the stau can be long-lived if the mass difference with the neutralino is smaller 
than Am < 2 GeV. Interestingly, the neutralino-stau coannihilation region, where the 
WMAP relic density can be obtained, also occurs for small mass differences. We have 
studied this possibility within the CMSSM, finding that once more the largest results 
for the stau flux are obtained for small values of the gaugino mass parameter, i.e., 
for light staus, and increase when the rest of the spectrum (in particular the squark 
masses) is also light. This favours low values of tan/5. For tan/3 = 10, and using the 
MPR flux, an upper bound of approximately 2 events yr _1 is found for staus with a 
mass of ~ 100 GeV, which decreases by a factor 2 for tan/5 = 50. When the 
WMAP constraint is imposed on the neutralino relic abundance, compatibility with 
viable neutralino dark matter reduces the allowed parameter space to a small range 
of stau masses. For example, with tan/3 = 10 one finds ~ 300 GeV, and the 
predicted stau flux is 0.1 events yr^ 1 . When tan (3 increases the region compatible is 
shifted to heavier staus and the resulting flux decreases significantly. Overall we have 
observed that the stau flux obtained in these supergravity scenarios is generally smaller 
than those of low energy supersymmetric analyses. This implies that, although stau 
events are distinguishable from the dimuon backround, in the most optimistic scenarios 
several years of data from IceCube are necessary to claim a positive signal. 
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Finally, we have investigated the case of a general supergravity theory in which 
the structure of the soft parameters is non- universal. We have observed that certain 
choices of non-universalities which lead to wino-like charginos and a decrease of the 
squark and slepton masses can account for a moderate increase of the resulting stau 
flux. More specifically, through a decrease of both the gluino and wino mass parameters 
with respect to the bino mass and a decrease of the soft mass for the H u Higgs at the 
GUT scale we have shown that the stau flux can increase by approximately a factor 
three with respect to the results in the CMSSM. 
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A Sfermion production from neutrino inelastic scat- 
tering 



The contribution from chargino exchange along a t-channel comprises the diagrams 
shown in the upper row in Fig. [7J They lead to the following cross-sections, 
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(A.12) 



with notation for Z matrices in agreement with that of [IE], where Z + and Z_ are 2x2 
chargino mixing matrices, with k — 1, 2, and Zjj, Zd and Z^ are 6x6 squark mixing 
matrices. Capital indices I and J are family indices running from 1 to 3, whereas 
indices % and j stand for mass eigenstates and run from 1 to 6. Summation over all 
indices is implied. 

For the diagrams with neutralino exchange along a t-channel, shown on the lower 
row of Fig. [7J we make the following definitions, 
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The total cross-section for interaction with quarks, expressed in terms of mass eigen- 
states, reads 
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(A.13) 

These parton-level cross-sections are then convoluted with their corresponding Parton 
Distribution Functions, which we extract from the Cteq6PDF package |36j. 

To calculate the fraction of the parent's energy carried by the stau, we use the 
general formula 
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Figure 7: Feynman diagrams describing the different neutrino inelastic interactions. 

where n is the number of intermediate states, i = corresponds to the parent squark 
or slepton, and i = n corresponds to the stau. 
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